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The two-component approach to the one-dimensional Dirac equation is apphed to the Woods-Saxon 
potential. The scattering and bound state solutions are derived and the conditions for a transmission 
resonance (when the transmission coefficient is unity) and supercriticality (when the particle bound state 
is at i^ = —m) are then derived. The square potential limit is discussed. The recent result that a 
finite-range symmetric potential barrier will have a transmission resonance of zero-momentum when the 
corresponding well supports a half-bound state at E — —m is demonstrated. 

INTRODUCTION 

o| 

^^ , There is a well-known theorem for low momentum scattering in the Schrodinger equation in one-dimension by an 

^^ ■ even potential well B, H: as momentum k tends to zero, the reflection coefficient L{k) tends to unity unless the 

potential V{x) supports a zero energy resonance 0. In this case L{k) — > and correspondingly the transmission 

coefficient T{k) — * 1. Bohm B calls this a transmission resonance. Recently we p] have generalised this result to the 

Dirac equation. Since the Dirac equation covers anti-particle scattering as well as particle scattering, the generalisation 

p^ . gives two distinct results since the fc ^ limit in the Dirac equation corresponds both to particle states where the 

energy E — m and anti-particle states where E — —m and m is the particle mass. The result that anti-particles 

can have transmission resonances when they scatter off potential wells is equivalent to particles having transmission 

^ ■ resonances when scattering off potential barriers. This is itself related to the result on barrier penetration found by 

^^ Klein |^ and now called the Klein Paradox. 

Our result B shows that transmission resonances at fc = in the Dirac equation occur for a potential barrier 

V = Uc{x) when the corresponding potential well V = —Uc{x) just supports a bound state a.t E = —m: this is called 

(^ . a supercritical state. While transmission resonances for a square barrier in the Dirac equation have been known for 

some time M] and their relationship to supercritical states for a square well was pointed out more recently S we are 

f^ , not aware of any other analytic solution of the Dirac equation for which they can be demonstrated. Nevertheless Dosch 

et al g showed in 1971 that a transmission resonance does exist for the Woods-Saxon potential, which is a smoothed 

out form of the square well/barrier. In this note, we take this example further. We solve the Dirac equation for the 

Q^' Woods-Saxon potential well to demonstrate the supercritical states and we find complete solutions for the reflection 

qj ' and transmission amplitudes for scattering off a Woods-Saxon potential barrier. We then are able to demonstrate the 

,^ . correspondence between the supercritical states and the transmission resonances analytically. We also consider the 

^ limit where the Woods-Saxon potential becomes a square well/barrier. 

H THE ONE PARTICLE DIRAC EQUATION IN ONE DIMENSION 

In the one-dimensional Dirac equation, solutions can be greatly simplified by adopting a two-component approach; 
both the positive and negative energy solution states arc retained without the added complication of spin. Starting 
with the rclativistic free particle Dirac equation (h — c — 1): 

i^r-^-m)i, = (1) 

In the presence of an external potential V{x) and taking the gamma matrices jx and 70 to be the Pauli matrices iax 
and az respectively, the Dirac equation in one dimension can be written as 

{<7x—-{E~V{x))<J,+ml)i;{x)^0 (2) 

ax 

The four-spinor, ip, is decomposed into two spinors, ui and U2, so that 






:a 



Thus the problem is to solve the coupled differential equations: 

u'l = -(m + E - V{x))u2{x) 
u'2 = -(m - E + V{x))ui{x) 

Following a similar procedure to that used by Fliigge [ |lO| , introduce the following combinations 

4i{x) — ui{x) + iu2{x) , x(2;) = ""1(0;) - iu2(a;) 

Substituting these into (4) and re-arranging gives: 

(/)'(x) — —imx{x) + i{E — V{x))(l){x) 
x'ix) — im(j){x) — i{E ~ V{x))x{x) 

The two components, (/"(a;) and x{x)j satisfy: 

0"(x) + [{E ~ V{x)f -m^+ iV {x)]4>{x) = 
x"{x) + [{E - V{x)f - w? - iY'(x)\x(x) = 

In the following the full solutions for 0(x) will be presented. In order to establish x(a:), use will be made of (6). 



(4) 
(5) 



(6) 

(7) 



(8) 
(9) 



THE WOODS-SAXON POTENTIAL 



The Woods-Saxon potential is defined as 



V{x) = W 



e{-x) 



0{x) 



1 + Q-o-i'-^+L) \ _|_ ga(a:-L) 



(10) 



with W real and positive for a barier or negative for a well; a and L are real and positive. 0{x) is the Heaviside step 
function. 
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FIG. 1. The Woods-Saxon potential barrier for L = 10 with a = 5 (solid line) and a — 0.5 (dotted line). 



At this stage it is worth mentioning that we will be interested in potentials where aL ^ 1; from the above figure 
it can be seen that, for this condition, the potential has a less pronounced cusp at a; = - the potential now closely 
resembles a square barrier with smooth walls. As stated in Q this docs not introduce any essential physical restriction 
on the problem and is significant only in that it allows exact solutions to be established (albeit in an approximation). 



SCATTERING STATES 



First, consider the scattering solutions for a; < with \E\ > m. On making the substitution y — —e °(^+^)^ 
equation (8) becomes 



2 "^ 

a y-r- 

dy 



dy 



E- 



W 



iayW 



= 



(11) 



Splitting off fitting powers of y and(l — y) by setting (f>L = J/^(1 — y) ^ f{y) and substituting into the above equation 
reduces it to the hypergeometrical equation 



y(i - v)ny) + [(1 + 2/i) - y(i + 2m - 2A)]/'(y) -[{ii-y- A)Oi + y- X)]f{y) = 

where the primes denote derivatives with respect to y and the following abbreviations have been used 

ip ik iW 

M=— , 1^ ^ — , A= 

a a a 



(12) 



p^ ^{E~ Wf ~ m^ 



E — m 



(13) 



Note that as we are considering scattering states, \E\ > m which ensures that k is real, and W is real and positive, p 
is real ior m < E < W — m (the Klein range) and E > W + m and imaginary for W — m<E<W + m. Of principal 
interest are the energies which lie in the Klein range and potentially lead to Klein tunnelling - where fermions can 
tunnel through strong potentials without exponential suppression [pj|. Equation (12) has the general solution 



f{y) = Di 2/-2M 2F1 {-fi - u - X, -n + ly - X,l - 2fi;y) + D2 2Fi{fi - ly - X, fi + ly - X,l + 2/i; y) (14) 



So 



My) = Di 2/-^(l - y)-^ 2Fl{-^l - I. - A, -/i + // - A, 1 - 2/i; y) 

+D2 y^{l - y)-\Fi{fi - J. - A, /i + J/ - A, 1 + 2/t; y) (15) 

For this to be a physically acceptable solution to the problem, it must satisfy the appropriate boundary conditions as 
a; — > — 00. The solutions as x —^ —00 => y —> —00 can be determined using the following formula for the asymptotic 
behaviour of the hypergeometric function [n2| : 



, , , r(c)r(6 - a) , , „ r(c)r(a -b) , , . 

2F,ia,b,c;y) = -^ L{-y)- + -^ M-y)-' 

1 (6)1 (c — a) 1 (a)l (c — 0) 

and noting that in the limit x -^ —00, (—y)^'^ -^ ^±ik{x+L) ^ Therefore in this limit, 4)l{x) can be written 

lim <J3l{x) = Ae''^'(="+^) + Be-'''^''+^^ 



(16) 



(17) 



From equation (6) the other component, xi^) is 

X{x) = — 



[i{E - V{x))<P{x) - <P'{x)] 



Substituting equation (17) into the above gives us 

lim xl{x) = A i I e 



where in both cases 

A = Di 
and 



r(i - 2/i)r(-2z^) 



T{-^i - ly - X)T{1 - fi - ly + X) 



ik(x+L) I r> 



D2 



E + k 



-ik(x-\-L) 



r(i + 2/i)r(-2z/) 



T{p - ly - X)T{1 + fi - ly + X) 



(18) 



(19) 



(20) 



B = D r(l - 2^,)T{2y) ^_„,^ ^ r(l + 2^i}T{2v) 

The choice of combinations of the wave function components (5) can be re-written : 

ui{x) = \{(l){x)+x{x)) , U2{x) ^ h^{x) ~ x{x)) (22) 

2 2,1 

Upon substitution of equations (17) and (19) into the above it can be seen that the wave function, tpix), comprises 
of an incident and reflected wave far to the left of the barrier which is the desired form to establish reflection and 
transmission amplitudes. 

Now consider the solutions for x > 0. The analysis will differ slightly from ||9| by making a more appropriate 
substitution to lead to the desired transmitted wave function far to the right of the barrier. This substitution will also 
lead to the correct wave functions when the bound state solutions are considered. On choosing z~^ = 1 + e"'^^^^^ 
equation (8) becomes 



a z(l — z)-— 
dz 



z(l-z) 



+ [iE-Wzf -m^ -iaz{l-z)W]<j>R^O (23) 



dz 

Putting (j)fj = z^^{\ — z)^'^g{z) and substituting into the above gives the hypergeometrical equation 

z(l - z)g"{z) + [(1 - 21,) ~ z(2 - 2m - 2i/)]g'(z) - [(1 - ^ - j. - X)(-^ - ^ + X)]g(z) = (24) 

The general solution to the above is 

g{z) ^ di2Fi{l- fi- ly - X,-^- iy + X,l-2iy;z) + d2 z'^" 2-Fi(l - ^i + i, - X, -^i + v + X,l + 2i'; z) (25) 

So 

(l)R = diz"''(l - z)~'^2Fi{l - i^L - V - X,-^. - V + X,l - 2v; z) + 

da z^il - z)-^ 2Fi{l- ^i + v-\-^l + v + \l + '^v\ z) (26) 

Also as X ^ oo, z ^ and z~'^ —> qI-H^-l) ^ Therefore in order to have a plane wave travelling to the right as a: ^ oo, 
da = 0. So 

0^ = di z-'^(l - zy 2Fi{l - fi - v ~ X, -fi - ly + X,l - 2iy; z) (27) 

lim ^R = di e-^'=^e^'=" (28) 



and 



whilst the other component 

lim XR = di (^^) e-'^'^e'"'-' (29) 

X— "-oQ \ m J 

In order to find the energy eigenvalues, the wave functions ^^ and (J)r must be matched ata; = 0. Asa;— >0,2/-^-0 
and z ^ 1 {aL ::^ 1), so 

and 

,, ^ ., ,-..,..,. r(2rir(i 2.) ^ ^_ e-..e-^.. ,^ r(-2rir(i 2.) 

1 (1 + /i — z/ — A)l (/i — J/ + A) 1 (1 — /i — i^ — A)l (— /i — V + X) 

where use has been made of the following continuation identity for the hypergeometric function in (pR: 

„ / . N r(c)r(c — a — 6) ^ , , 

2F1 a, 6, c; z = ,-r., rr 2^^1(0, &, a + 6 - c + 1; 1 - z + 

1 (c — a)l (c — 0) 

(1 - ^r^'^^'^-^^^^^^^fp^ 2Fi(c - a, c - 6, c - a - 6 + 1; 1 - z) (32) 

T{a)T[h) 



Comparing the coefRcients of e^°'^^^ and eliminating di gives 



£l = ,-2».u 4aLM r(2M)r(i -li-y- A)r(-M - ^^ + A) 

The electrical current density for the one-dimensional equation (2) is defined as 

j = ilj{x)-fx'4^{x) = -ip(x)'^ a2'tp{x) 

— i {UiU2 ~ U2U1) 

= ^ (10(^)1' -Ix(2^)r) 



(33) 



(34) 



where use has been made of the choice of combinations for the wave functions (5). The current as a; — > — cxd is 
JL ~ jin — jrefi whcrc jin is the incident current and jrefi is the reflected current. Similarly as x — > cx3 we have that 
the current jr — jtrans where jtrans is the transmitted current. Substituting equations (17), (19), (28) and (29) into 
(34) we find that 



3L= \A\'^{E-k)-\B\''^{E + k) 
m= \di\'A{E~k) 



12 fc 
m 



(35) 



From the conservation of charge we have that jl = JR which together with the relection coefficient, R, and the 
transmission coefficient, T 



R 



Jrefi 
Jin 



\B\ 



E + k 



|v4|2 \E-k 



T = 



Jtrar, 



\A? 



we obtain the unitarity condition: 



R + T = 1 



(36) 



(37) 



For the purpose of comparison with M the following expressions for the reflection and transmission amplitudes are 
useful : 



m + E + k 



m + E -k , 
The reflection amplitude is found to be 

_ 6-2*'=-^ {m + E + k)B{2v,-n-v-X) 



-2ikL^ 
A 



-2ikL 



A 



n [m + E -k)B{-2v,l- ^i + v + \) 



:„i 5(2^, -^i-v + X)B{2n, -11 + v + A) 



^iipL 



where 



17 == 1 - e"*"-^^ 



[jx + uf-X 
in - pf - A2 



21 



B{-2fi, ^i + v - X)B{-2^, ^i-v -X) 



B^^{2^l,-^l-v + X) 
B2(-2^,/i-i^-A) 



and B{a^ b) = J^, , , n is the Beta function. Note that this equation yields identical results to the r given in 
transmission amplitude is found to be 



t 



n 



Afiv 



r2(-2^)r2(-2i/) 



(38) 
(39) 

(40) 
The 

(41) 



Using equations (39) and (41), the unitary condition (37) can be established; the algebra is quite involved and some 
care must be taken for the two cases where p is real and imaginary. 




W-m 




E-m 



E+m 



W 



FIG. 2. The transmission coefficient for the relativistic Woods-Saxon potential barrier. Tlie left hand plot illustrates T for 
varying energy, E, with L = 10, a = 5, m = 0.4 and W = 1.2. The right hand plot is for varying barrier height, W , with 
L = IQ, a = f>, m = 0.2 and E — 2m. Both plots illustrate the tunnelling without reffection predicted by Dosch, Jensen and 
Miiller. 



In order to establish the condition for a transmission resonance, T = 1, to occur, it is more convenient to look at 
the reflection amplitude for r = 0. In this instance, only the square-bracketed term on the right of (41) needs to be 
considered; although poles of the gamma-functions could make r zero, /x, v and A are all purely imaginary for energies 
in the Klein range and do not give the negative integer required for this to happen. Consequently, transmission 
resonances occur when 



^ _ ^4^pL B{2ii, -ii-u + A)g(2/i, -ii + u + X) ^ ^ 



(42) 



BOUND STATES 



In order to study the bound states, use can be made of the wave function derived for a; > 0, but the analysis can 
be simplified for a; < by making the substitution j/~^ = 1 + 6^"^^+^^ which leads to the following equation: 



«^J/(l-y)^ 






+ [{E + Wyf -m^ - iay{\ - y)W\ 







(43) 



where W -^ —W in equation (10) for potential wells. Putting 0^ = 2/'^(l — yVh{y) leads to the hypergeometric 
equation 



2/(1 - y)h"{y) + [(1 + 2cr) - j/(2 + 2cr + 2-i)]h'{y) - [(1 + a + 7 - A)(cr + 7 + \)]h{y) = 
where 

a = - , K^m^-E^ , 7=^ , p'^^iE + Wf-m" 
a a 

Asa;^— oo,y^Oso choose the solution 

h{y) = 2i^i (1 + (T + 7 - A, o- + 7 + A, 1 + 2cr; y) 
and therefore 

4)L = A'y%l - yy2Fi{l + a + 7 - A, a + 7 + A, 1 + 2a; y) 



(44) 



(45) 



(46) 



(47) 



For X > 0, use equation (26) with W -^ —W so A ^ —A and /i ^ 7 and also k ^> ik so v ^> —a then choose the 
solution 



B'z^il - z)-^2Fi{l + a - 7 + A, cr - 7 - A, 1 + 2a; ;j) 



(48) 



Once again, in order to find the energy eigenvlaues these two wave functions must be matched at a; == where both 
y, 2: — > 1 (aL ^ 1). By making use of the continuation formula (32) the component (j) can be written 



and 



j^l^-iax^-iaL 



<pR -^ B'e-^''^e'^"-^- 



r(-27)r(l + 2a) 



r(l - 7 + tJ - A)r(-7 + o- + A) 
r(27)r(l + 2a) 



+ A'e'^'^^e''"'^ 



r(27)r(l + 2a) 



5'e7aXg-7aL_ 



r(l+7 + a-A)r(7 + a + A) 
r(-27)r(l + 2a) 



r(l + 7 + a + A)r(7 + a - A) r(l - 7 + a - A)r(-7 + a - A) 

Comparing terms in e '''"^ and chminating A' and B' ultimately leads to 



g(-27,7 + a-A)^ 
B(27,-7 + a + A)2 



(a + 7)2 - A2 



So 



g(-27,7 + a-A) _^^2-,aL 
B(27,-7 + a + A) 



= ±e" 



'(a-7)2-A2 
(a + 7)2 - A2 



(49) 



(50) 



(51) 



(52) 



where the even solutions are determined by the positive square root and the odd solutions by the negative square 
root. These equations need to be solved numerically to find the energy eigenvalues for the bound states. It is also 
possible to plot the upper and lower components of the bound state wave function tfj{x): 




FIG. 3. The upper (left) and lower (right) components of the first even bound state wave function for the sub-critical 
Woods-Saxon potential well where a = 10, L = 2, m = 1 and W = 2. The energy of the bound state is -B = 0.759. 



SUPERCRITICALITY 

As the potential well deepens for increasing W , the energy eigenvalue of any given bound state will also decrease. 
When this energy reaches E = —m, the bound state merges with the negative energy continuum and the potential is 
said to be supercritical. When E —> —m, k —^ and consequently a — > 0. Writing 



7c = -Pc 
a 



the energy eigenvalue equation (51) becomes 

S(-27c,7c 



Pc 



A)^ 



S(27„ -7, + A)^ 



14^2 _ 2mW 



= g47cai 



(53) 



(54) 



Once again the even and odd supercritical energy eigenvalues can be determined by taking the positive and negative 
square roots respectively. 



SQUARE WELL LIMIT : BOUND STATES 

Using 



Bia,b)^{^-p)(^-±^^)Bia + l,b+l) (55) 



we find that 



B(T27, ±7 + - T A) = T (lI±fiMii2±^i^5(i ^ 2,, 1 ± , + . T A) (56) 



Therefore in the square weli hniit a— >oo,7,cr, A^O and 

B(-27,7 + g- A) (-7 + g- A)(-7 + g- + A) 

5(27, -7 + cr + A) ^ (7 + 0-- A)(7 + cr + A) 

where we have used B{1, 1) = 1. Substituting this into (51) we obtain 



e 



4ip' L 



(g - 7)^ - A^ _ (k - ^y)2 + W^ 
(<T + 7)2 - A2 ^ (k + ip')"^ + W^ 



(57) 



(58) 



Rationahsing we find that 



2mW mW 

Choosing the positive root and soiving for the real and imaginary parts one eventually obtains 



^-"- C''/!.:"" -^! (58) 



2Kn' 
^«" 2p'i = ^5 I---J (60) 



which after much laborious algebra gives for the even solutions 



, mW + K^-EW (m-E){E + W + m) 

tanp L = ^ \ -, Z-, r (61) 

^ Kp' \ {m + E){E + W -m) ^ ' 



and for the odd solutions 



, mW-n'^+EW {m + E)(E + W-m) 

'"'''''' = 1^' = ~'^ {n.-E){E + W + m) ^^'^ 

(these are precisely the equations for even and odd bound states in the square well [g| with W — V). 

ZERO MOMENTUM RESONANCES AND SUPERCRITICALITY 

It was first pointed out in [|l3| that the conditions for supercriticality and zero-momentum resonances were the same 
for the square, Gaussian and Woods-Saxon potentials. Indeed for the square potential well, V{x), where V{x) = 
for |a;| > a and V{x) = —U < for |a;| < a, the condition for supercriticality is 2p'a = Nw where p'^ = C/^ — 2mU 
|0|. The Dirac equation (4) is invariant under charge conjugation: that is to say under the transformation 

E -^ —E U -^ — [/ ui — > M2 U2 — > ui (63) 

Consequently the condition for a supercritical particle a,t E = —m in a square well is the same as that for a supercritical 
antiparticle at -E = m in a square barrier. It can also be shown that for a transmission resonance to occur, 2pa — Nt: 
where p^ ~ {E — V)'^ — rn^ = k"^ + V'^ — 2VE. In the zero-momentum limit, fc ^ 0, this is seen to be identical to 
the condition for a supercritical antiparticle. In other words when a potential well of finite range is strong enough 
to contain a supercritical state, then a particle of arbitrarily small momentum will be able to tunnel right through 



the potential barrier created by inverting the well without reflection [pi. From (42), the condition for a transmission 
resonance to occur for the Woods-Saxon potential is 

When E ^ m (for a low- momentum particle), we find that p^ — > W^ — 2mW — Pc so that jjl ^> ^c and also 
fc ^ ^ i^ ^ 0. So (64) becomes 



PsL _ 



B{-2jc,lc~Xy 
B(27„-7e + A)2 



(65) 



As expected, this is precisely the condition required for the potential barrier to be supercritical (54) once the 'flip- 
ping' procedure on the potential well as described above is considered. The components of the zero-momentum 
resonance/half-bound state wave function, ^{x)^ can be plotted and have the following appearance: 




FIG. 4. The upper (left) and lower (right) components of the first zero momentum wave function for the supercritical 
Woods-Saxon potential barrier where a = 10, L = 2, m = 1 and W — 2.274. 

Comparison with the zero momentum wave functions for the square and Gaussian barriers H highlights the existence 
of two turning points which occur in the top component of the wave function for the smooth Gaussian potential only. 
Once again these correspond to the points ixjf where V{±xk) ^ E + rn = 2m at zero momentum and are not 
manifest for the square barrier whose walls are discontinuous. 
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